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ABSTRACT 



^ ' Aims. We quantify the effect of gravitational redshift on emission lines to explore the transition region from the Newtonian to 
the Einsteinian regime. With the emitting region closer to the Kerr black hole, lines are successively subjected to a stronger 
gravitationally induced shift and distortion. Simulated lines are compared to broad, optical emission lines observed in Mrk 110. 
Methods. We simulate relativistic emission line profiles by using Kerr ray tracing techniques. Emitting regions are assumed to 
be thin equatorial rings in stationary Keplerian rotation. The emission lines are characterised by a generalized Doppler factor 
or redshift associated with the line core. 
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, Results. With decreasing distance from the black hole, the gravitational redshift starts to smoothly deviate from the Newtonian 
Doppler factor: Shifts of the line cores reveal an effect at levels of 0.0015 to 60 % at gravitational radii ranging from 10 J to 2. 
This corresponds to fully relativistic Doppler factors of 0.999985 to 0.4048. The intrinsic line shape distortion by strong gravity 
i.e. very asymmetric lines occur at radii smaller than roughly ten gravitational radii. 
Q ' Conclusions. Due to the asymptotical flatness of black hole space-time, GR effects are ubiquitous and their onset can be tested 
S— i observationally with sufficient spectral resolution. With a resolving power of ~ 100000, yielding a resolution of ~ 0.1 A for optical 
C/} and near-infrared broad emission lines like H/3, Hell and Paa, the gravitational redshift can be probed out to approximately 
75000 gravitational radii. In general, gravitational redshift is an important indicator of black hole mass and disk inclination as 
recently demonstrated by observations of optical lines in Mrk 110. Comparing our simulated lines with this observations, we 
• ^ ■ independently confirm an inclination angle of 30 degrees for the accretion disk. Redshift deviations induced by black hole spin 
p\i ' can be probed only very close to the black hole e.g. with X-ray iron lines. 
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1. Introduction 

Active galactic nuclei (AGN) such as Seyfert galaxies 
and quasars are powered by accreting supermassive black 
holes (SMBHs) following the standa rd model that has 
been developed over fou r decades ijLvnden-Belll [l969; 
iLvnden-Bell fc Reed Il97l|) . The m asses of SMBH s lie in 
the range 10 6 to 10 10 M Q , see e.g lNetzerl l)2003l) . In the 
standard model, clouds moving in the gravitational poten- 
tial of the black hole are photoionized by the central AGN 
continuum, thereby producing Doppler broaden ed emis- 
sion l ines with widths of typically 10 3 -10 4 kms -1 ijWoltierl 
Il959l as pioneering studies). The region where the broad 
lines originate is usually referred to as the broad-line re- 
gion (BLR). The scale of the BLR is believed to be 10 15 
to 10 17 cm, corresponding to ~ 10 3 to ~ 10 5 r g or 0.6 
to 60 light days for a 10 7 M Q black hole. Here the grav- 
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itational radius is defined as r g = GM/c 2 with Newton's 
co nstant G, vacuum spee d of light c and black hole mass 
M. lRobinson et all ( 199ffl have presented complex models 
involving spherical or disk geometries for the BLR as well 
as rotational and radial cloud kinematics. They studied 
the influence of continuum variability on line profiles and 
found diverse line shapes exhibiting spikes, bumps and 
shoulders though in a non-relativistic regime. Additional 
velocity components in the BLR caused by disk winds 
l|Konigl fc Kartieil994|) and a radial component have been 
suggested from accretion theory and radia l velocity maps 
of the narrow-line region l|Ruiz et alJl200lft . However, the 
detailed structure and vel ocity field of the BLR remain 
unclear l|Collin et al J I2006) . 

The broad-line clouds respond to variations in the central 
photoionizing continuum as suggested by strong correla- 
tions between H/3 response times and non-stellar optica l 

H 12002ft . 



continuum fluxes, r ccnt oc vTuv IjPeterson et al.l 
This phenomenon is exploited in reverberation mapping 
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techniques to deter mine both the scale of the BLR and 
the black hole mass jglandford fc McKedll983 IPetersorJ 
llQQfllKaspi et alJl2000l) .~ 

The idea that gravitational redshift may influen ce opti- 
cal lin es causing line asymmetries was raised by iNetzerl 
JHZZ|). In the Seyfert-1 galaxy Akn 120, a slight redward 
displacement of the H/3 line was reported, amounting to 
Az ~ 0.0013, interpreted a s the result of gravitational red- 
shift l|Peterson et alJl985j) . However, such effects may also 
arise from attenuatio n of the BLR or light - travel time ef- 
fects, as discussed bv lPeterson et all l)l985|) . Similar stud- 
ies that assume that observed effects arc a result of gravi- 
tational redshift have been done for a quasar sample where 
the SMBH mass of QSO 0026+129 could be roug hly esti- 
mated to be 2 x 10 9 Mq l|Zheng fe Sulentid[l9901 ): this is 



still the current value within a factor of 2 



Czernv et al.1 



2004) . Recently, several BLR optical emission lines in the 
narr ow-line Seyfert-1 galaxy Mrk 110 were investigated 
(see lKollatschnvll2003l K03 hereafter). In that work, Ha, 
H/3, HeIA5876 and HeIIA4686 emission lines were found to 
possess a systematic shift to the red, with higher ioniza- 
tion lines showing larger shifts as expected in a BLR with 
stratified ionization structure. 

In this paper, we study the gravitational redshift over a 
large range of distances from the central black hole. We 
quantify the relativistic gravitational redshift on emis- 
sion lines until GR fades beyond the current observable 
limit. The investigation is carried out in a very general 
form by discussing the observed line profile as a func- 
tion of the generalized GR Doppler factor (g-factor) for 
Kerr black holes and an arbitrary veloc ity field of emit- 
ters, see e.g. lMuller fc Camenzindl ll2004 M04 hereafter). 
Pioneering work on relativistic spectra was performed by 
iGunninghaml (jl975ll using transfer functions. However, the 
considerations of the g-factors in this work were restricted 
to minimum and maximum values of g on infinitesimally 
narrow and thin stationary rings. Furthermore, the dis- 
tance range of interest for BLRs, 10 3 to 10 5 r g , has not 
been investiga ted in detail. 

iGorbinl ( 1997t) studied relativistic effects on emission lines 
from the BLR by assuming Keplerian orbits for the emit- 
ting clouds in a Schwarzschild geometry. It was found that 
line profiles decrease in both, width and redward centroid 
shift when the line emitting region moves away from the 
black hole. 

Our goal is to accurately quantify the effects of gravita- 
tional redshift in the vicinity of a Kerr black hole. After 
a very general consideration that holds for any classical 
black hole of arbitrary mass, a more specific treatment in- 
volving optical emission lines from BLRs is addressed. For 
the case study of Mrk 110, it is even demonstrated how 
the mass of the SMBH and the inclination of the inner 
disk can be determined. 



2. Method, analysis tools and model 

2.1. Relativistic ray tracing 

In contrast to Cunningham's work, emission lines are com- 
puted by ray tracing in the Kerr geometry of rotating 
black holes. Light rays emitted in the vicinity of the black 
hole travel to the observer on null geodesies in curved 
space-time, and in this work the observer is assumed to be 
located at r Q b s = 10 7 r g . The Kerr Black Hole Ray Tracer 
(KBHRT) maps emitting points in the equatorial plane 
of a Kerr black hole to points on the observer's screen. 
Spectral line fluxes are computed by numerical integra- 
tion over the solid angle subtended by the screen. All rel- 
ativistic effects such as gravitational redshift, beaming and 
lensing are included, but higher order images are not con- 
sidered. The complete solver has been presented in earlier 
work (M04). 



2.2. Analysing relativistic emission lines 

In the following, line fluxes are discussed as a function of 
the g-factor which is defined as 



— A cm /A, 



1 



obs 



l + z 



(1) 



where v and A denote frequency and wavelength, respec- 
tively, and the redshift is z. Emitter's and observer's frame 
of reference are indicated by subscripts 'em' and 'obs'. A g- 
factor of unity therefore corresponds to an unshiftcd line, 
whereas g < 1 indicates rcdshiftcd emission and g > 1 
blueshifted emission. Note that by using the g-factor, the 
simulated line profiles are discussed very generally without 
specifying a particular emission line. 

Analysis tools for relativistic emission lines as well as 
line classification schemes by morphology that were in- 
troduced by M04 (Sec. 8) are utilized here. A relativistic 
emission line exhibiting two Doppler peaks can be char- 
acterised by several quantities: 

— ffmin, minimum g-factor that defines the terminating 
energy at the red wing 

— g max , maximum <?-factor that defines the terminating 
energy at the blue wing 

— g rp , the g-factor associated with the red relic Doppler 
peak 

— <7bp) the g-factor associated with the blue beamed 
Doppler peak 

— DPS = gb p - 5rp, the Doppler peak spacing in energy 

— F max , maximum line flux 

— F rp , line flux associated with the red relic Doppler peak 

— Fb p , line flux associated with the blue beamed Doppler 
peak 

— DPR = Fb p /F rp , the flux ratio of both Doppler peaks 

The existence and strength of these parameters depend 
on the line shape. Low inclination angles of the emitting 
surface i.e. face-on situations with axial observers destroy 
the typical structure with two distinct Doppler peaks. 
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Based on relativistic emission line terminology, line 
morphologies can be classified as triangular, double- 
peaked, double horned, shoulder-like and bumpy shapes. 
Triangular and shoulder-like morphologies lack a red 
Doppler peak. Bumpy morphology even lacks a distinct 
blue beaming peak either because a very steep disk emis- 
sivity suppresses emission at large disk radii or because 
the line originates too close to th e black hole. We want 
to stress here that iRobinson et alJ l|l990|) found a similar 
terminology for non-relativistic lines but the classification 
scheme for relativistic lines in M04 was established inde- 
pendently. 

In order to be able to characterize any line profile in- 
dependently of its morphology, and in order quantify the 
shift of the resulting line centroid, we define a new quan- 
tity, Score : 



9c 



(2) 



The g core parameter is thus associated with the line core 
energy, i.e. the energy associated with the flux weight of 
the whole line. As can be seen, it is evaluated numerically 
by multiplying each energy bin, g- u with the corresponding 
flux in that bin, F\, and summing over the line profile. 

Gravitational redshift in the weak field regime estab- 
lishes pure shifts of spectral features without changing 
their intrinsical shape. However, gravitational redshift in 
the strong field regime - strong gravity - produces re- 
markable distortions of spectral shapes if the rest frame 
feature is compared to its analogue in the observer's frame. 
Distortion is a key feature of relativistic sp ectra exhibit- 
ing very skewed and asymmetric line profiles l|Fabian et alJ 
Il989t IPopovic et alJIlQQfil iTanaka et alJll99fih . These ef- 
fects are important only very close to the black hole. 
Here we investigate how gravitational redshift changes the 
mode i.e. Einsteinian gravity transmutes to Newtonian 
gravity for emission regions moved away from the black 
hole. 

We also define the half-energy radius as the radius 
where a given g-factor is reduced to exactly 1/2 its original 
value. More precisely, Rj denotes the half-energy radius 
associated with j e {g min , ffmax, Score, ffrp, 3b P }- Due to 
Eq. JU the observed energy of the radiation associated 
with the specific g is exactly one half of the emitted energy 
in the rest frame. As a measure of strong gravity i? grp is 
chosen as the radius where the g-factor associated with 
the red Doppler peak is exactly 1/2. This is a suitable 
choice because strong gravity deforms the red line wing in 
an extraordinary manner. 

2.3. Emitter model: Rendering parameters and 
emissivity 

The aim of this paper is to study gravitational redshift 
effects. Therefore, it is of particular interest to avoid 
blueshift effects that would blur or cancel the redshift. An 
easy way to switch off blueward displacements is tilting 
the emitting region to a face-on situation: axial observers 



with inclination angle i = 0° to the emitting area have 
no relative motion along the line of sight to the emitter. 
In this case, the radiation is only affected by gravitational 
redshift because no motion is directed out of the plane. 

For numerical reasons an inclination angle of i = 1° 
is chosen which is close enough to the face-on situa- 
tion. Other parameters are: Kerr black h ole rotating at 
Thome's limit a/M = 0.998 i|Thornell974h and a prograde 
Keplerian velocity field of emitters, J1k = y/M/(Vr^ + 
ay/~M). We consider here stationary thin rings with emis- 
sion peaking at i? pca k- This is established by rendering a 
disk and shifting a Gaussian radial emissivity profile over 
the disk 



e(r) cx exp 



(r - -Rpeak)' 



(3) 



as introduced by M04. The parameter er r controls the 
width of the Gaussian or the size of the emitting region 
and is chosen to be 0.2. The peak radius, i? pca k, is given in 
units of gravitational radius. A localized Gaussian emissiv- 
ity set in this way mimics a thin and narrow luminous ring 
with ~ 1 r g distance between inner and outer edge. The 
Gaussian guarantees a smooth but steep decay of emission 
at the ring edges. Note that using r g as a natural scale al- 
lows for applications to any classical black hole including 
also stellar and intermediate-mass black holes. 

The technique is now straightforward: The ring is 
shifted from large distances in asymptotically flat space- 
time in the direction towards the black hole where space- 
time curvature becomes strong. In the present work, we 
assume that the radial range of interest is 2 r g < i? pca k < 
100000 r g . For each simulated emission ring we determine 
the line core energy g corc by applying Eq. and comput- 
ing the associated line redshift, z core , from Eq. Q. These 
quantities allow for quantifying the gravitational redshift 
as a function of distance to the rotating black hole. 

Further, we use the ring emitter model as a simple 
model for the broad-line region in AGN. In this case, 
the BLR clouds are distributed in an equatorial plane 
and follow Kepleria n orbits around the black hole as 
proposed els ewhere ijNetzerl Il977t IRobinson et all ^990; 

Depending on radial distance to the cen- 
ter and orientation angle to observer, the BLR emission 
is influenced by Doppler effect, gravitational redshift and 
beaming to different extent. We show in Sec. 01 that it is 
possible to fit optical data to this simple, flat Keplerian 
BLR model and that plausible inclination angles of the 
inner disk can be deduced. 

3. Gravitationally redshifted emission lines 

3.1. Redshifted line cores 

As outlined in the previous section the line cores are com- 
puted by ray tracing and their dependence on radial dis- 
tance to the black hole is analysed. Fig.^displays the core 
redshifts of emission lines, z cmo , as a function of i? pca k for 
rings inclined to i = 1°. The redshifted lines are dominated 



— 0.001 




A. Miiller and M. Wold: Signatures of Gravitational Redshift 

Table 1. Radial redshift dependence for i = l c 



cr 



Fig. 1. Radial dependence of redshift z of line cores (filled 
circles, left y-axis) and strength of gravitational redshift 
effect (triangles, right y-axis). Inclination angle amounts 
1°. 



by gravitational redshift by construction. The filled circles 
show that the redshift approaches z — > 0, i.e. that g — ► 1, 
at distances of a few thousand gravitational radii from the 
black hole. This is the regime of nearly flat space-time 
and Newtonian physics. But approaching the black hole, 
space-time curvature becomes more significant: z grows 
rapidly and g approaches zero. 

The triangles illustrate the strength of the gravita- 
tional redshift effect in an alternative way. The scale of the 
axis on the right-hand side is computed as (1— g CO rc) x 100 
so that a value of g = 1 corresponds to 0% effect and 
g = to a 100% effect, i.e. the gravitational redshift at 
the event horizon of the black hole. The error bars result 
from uncertainties in determining g-factors. Ray tracing 
simulations with different numerical resolutions in both 
disk resolution for rendering and spectral resolution for 
line computation yield slightly deviating values for <? core 
(and other quantities in general). 

Gravitational redshift can alternatively be visualised 
by plotting the line core energies, g C ore, as a function of 
the peak radius for each ring. Fig. |3] shows the result for 
i = 1° and i = 75°. Highly inclined rings exhibit strong 
blueshift effects overlapping the redshift. As a consequence 
the core ^-factors of the i — 75° dataset never drop below 
0.8 for radii larger than the marginally stable orbit. So, 
comparison of both orientations demonstrates that highly 
inclined rings are not appropriate to study gravitational 
redshift as pure shifting effect. It is shown later that high 
inclinations are well-suited for probing strong gravity. 

In Tab.n we show the results of the Kerr ray tracing 
simulations for nearly face-on rings, i = 1°, in terms of 
numerical values for <? coro , z CO rc and GR effect as a func- 
tion of distance to a Kerr black hole with a/M = 0.998. In 
principle, Tab.^ilhistrates the die out of GR with increas- 
ing radius. However, it is important to note that according 
to the asymptotical flatness of GR black holes solutions, 



-^pcak [^g] 


^corc 


■2-corc 


GR effect [%] 


2 


0.404763 


1.470584 


59.52 


4 


0.649680 


0.539220 


35.03 


6 


0.755815 


0.323075 


24.42 


8 


0.813688 


0.228973 


18.63 


10 


0.849776 


0.176780 


15.02 


20 


0.924283 


0.081919 


7.57 


30 


0.949551 


0.053130 


5.04 


40 


0.962198 


0.039287 


3.78 


50 


0.969071 


0.031916 


3.09 


60 


0.975612 


0.024997 


2.44 


70 


0.978597 


0.021871 


2.14 


80 


0.980951 


0.019419 


1.90 


90 


0.982324 


0.017994 


1.77 


100 


0.984930 


0.015301 


1.51 


200 


0.992485 


0.007572 


0.75 


300 


0.994996 


0.005029 


0.50 


400 


0.996250 


0.003764 


0.38 


500 


0.997002 


0.003007 


0.30 


600 


0.997503 


0.002504 


0.25 


700 


0.997860 


0.002144 


0.21 


800 


0.998128 


0.001875 


0.19 


900 


0.998337 


0.001666 


0.17 


1000 


0.998503 


0.001499 


0.15 


2000 


0.999253 


0.000747 


0.07 


3000 


0.999502 


0.000498 


0.05 


4000 


0.999621 


0.000380 


0.04 


5000 


0.999696 


0.000304 


0.0304 


10000 


0.999846 


0.000154 


0.0154 


20000 


0.999923 


0.000077 


0.0077 


30000 


0.999949 


0.000051 


0.0051 


40000 


0.999961 


0.000039 


0.0039 


50000 


0.999970 


0.000030 


0.0030 


60000 


0.999975 


0.000025 


0.0025 


70000 


0.999978 


0.000022 


0.0022 


80000 


0.999981 


0.000019 


0.0019 


90000 


0.999983 


0.000017 


0.0017 


100000 


0.999985 


0.000015 


0.0015 



space-time curvature approaches zero only in the limit 
r — ► oo i.e. there is no finite distance at which the gravi- 
tational redshift vanishes exactly. Hence, Tab. could be 
generally continued ad infinitum. However, observability 
poses a limit to what is practical since the resolving power 
of a spectrograph, in the ideal case, constrains the amount 
of gravitational redshift (parameterized by z core here) that 
can be detected. Assuming a spectral resolution of 0.1 A 
for H/3 as is obtainable by instruments like UVES and 
CRIRES on the VLT, the corresponding critical value of 
the g-factor is g — 0.999979. In terms of velocity shift, this 
is w 10 kms -1 . As seen in Tab.^ this shift occurs at a ra- 
dius of ~ 75000 r s . Hence, we do not consider radii above 
100000 r g in Tab.Q] For supermassive black holes of 10 7 - 
10 8 M Q this radius corresponds to 0.05-0.5 pc, whereas 
for stellar-mass black holes of ~ 10 M Q it is 0.01 AU. 
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Fig. 2. Shift of line core energy in units of g-factor with 
distance to the black hole. Lowly inclined rings, i = 1°, 
are compared to highly inclined rings, i — 75°. Note the 
blueshift g > 1 in the latter case. 



Applying Eq. an optical Hell emission line with 
Aom = 4686 A in the emitter frame is gravitation- 
ally redshifted to 4686.1, 4686.7, 4693.0, 4757.7 A at 
100000, 10000, 1000, 100 r g . 

3.2. Line distortion by strong gravity 

In the previous section the gravitational redshift was in- 
vestigated as an effect that shifts the whole line core. We 
now focus on the behaviour of the intrinsic shape of rel- 
ativistic emission lines as a function of distance to the 
black hole. This is dubbed strong gravity as anticipated 
in Sec. 12.21 In the following, we tilt the emitting ring to an 
inclination of i = 75° in order to be able to study the char- 
acteristic broad profile with the two relic Doppler peaks. 
Fig. |3 shows how the line profile changes by tilting from 
i = 1° to i = 75°. 

Fig. 0| shows the effect of gravitational redshift on the 
i = 75° profile as a function of distance to the black hole. 
The line profile broadens and the line flux gets more and 
more suppressed as the emitting ring is moved closer to 
the black hole 1 . Eventually, the line decays and disappears 
at the event horizon. The red relic Doppler peak is shifted 
to lower energies as the ring approaches the hole, illus- 
trated in Fig. El The red peak flux is also more and more 
suppressed as can be seen in Fig.@] Close to the black hole 
the distortions are so strong that the red Doppler peak be- 
comes highly blurred and effectively vanishes in the line 
profile. At the event horizon the profile dies out and be- 
comes unobservable. Fig. \5\ can be used to read the half- 
energy radius associated with the red relic Doppler peak 



as defined in Sec. 12.21 At i = 75° this value is Rg T p — 5 r 



Fig. 3. Two example relativistic lines emitted from a 
highly inclined ring, i = 75°, and a face-on ring, i = 1°, 
both with the emission peaking at the same i? pca k = 
100.0 r g around a Kerr black hole with a/M = 0.998. 
The arrows indicate the core g-factors, g coro ~ 1.008 and 
g corc ~ 0.985 respectively for each particular case. 
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1 However, it must be noted that the line flux also gets 
suppressed because the effective emitting area decreases for 
smaller rings that have constant r out — n n in any case. 



g-factor 

Fig. 4. Line distortion by strong gravity. The emission 
of narrow rings peaks at the radii as denoted in the leg- 
end. There is no way to confuse the lines because line flux 
successively decreases with decreasing radius of maximum 
emission i? pca k- All rings follow Keplerian rotation around 
a Kerr black hole with a/M — 0.998 and all are inclined 
to an inclination angle of i = 75°. 



and documents that strong gravity is important only very 
close to the black hole. The distortion of the intrinsic line 
profile by gravitational redshift can be seen from g rp alone 
or by using another line criterion. A further line character- 
istic is the Doppler peak spacing (DPS) i.e. the energetic 
distance of both Doppler peaks (if available). DPS can be 
measured in units of energy, frequency, wavelength or - as 
has been done here for generality - in units of g. Fig. [|J] 
shows that this quantity does not remain constant as the 
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Fig. 5. Distortion of the red relic Doppler peak for rings 
satisfying i — 75°: With decreasing radius the red peak 
can be found at lower peak energies due to strong gravita- 
tional redshift. Very close to the black hole - r < 3 r g at 
this specific inclination - these distortions are such strong 
that the red Doppler peak is highly blurred and vanishes 
effectively in the line profile. 



Fig. 6. Energetic distance of red and blue Doppler peak 
at i = 75°. This Doppler peak spacing (DPS) is measured 
in units of g. Far away from the black hole both peaks 
approach significantly and the double-peaked line profile 
becomes very narrow: At ~ 3000 r g the peak difference in 
g only amounts to ~ 0.03. 



black hole is approached. DPS rises quickly so that the line 
is stretched. This phenomenon is directly but only quali- 
tatively visible in Fig. 0] Gravitational redshift causes an 
additional suppression in flux so that close to the black 
hole the emission line from an intermediately to highly 
inclined ring is asymmetric and skewed. 

We close this section with a comment on black hole 
detectability: Redshifted line cores as presented in Sec. 
13. II leave enough room for other gravitational sources than 
black holes; in contrast, spectral lines distorted by strong 
gravity in connection with a measured high compact mass 
support black hole candidates. 

4. Gravitationally redshifted optical emission lines 

4.1. Observations of NLS-1 Mrk 110 

In a recent work by iKollatschnvl l)2003|L broad optical 
emission lines from Ha, H/3, HeIA5876 and HeIIA4686 were 
found to display significant and systematic redshifts. By 
reverberation mapping, the distances of the emitting re- 
gions from the central continuum source were determined 
and a stratified ionization structure seen with Hell aris- 
ing closest to the black hole at a distance of ~ 490 r g . 
The observed shift of the Hell line was measured to be 
Az ~ 0.002, corresponding to g ~ 0.998. 

4.2. Inner disk inclination of Mrk 110 from Kerr ray 
tracing 

In this section, we follow the assumption that the observed 
optical lines of Mrk 110 are subject of gravitational red- 
shift and Doppler shifts and compare them with Kerr ray 



tracing techniques. It is aimed to determine parameters of 
the black hole-BLR system. The flat BLR model assum- 
ing line emitting rings as outlined in Sec. 12.31 is applied. 
Line redshifts are computed with the fully generalized GR 
Doppler factors g that were presented by M04 (see Eq. 13 
therein) . We note that in order to capture the correct g C ore 
value, a fine line binning is needed and hence a very high 
numerical spectral resolutions is used, Ag = 0.00001. 

In contrast to our redshift analysis we now have to 
allow for arbitrary inclination angles for the BLR in 
Mrk 110. A parameter space with inclination between 1° 
and 40° and radii in the range 100 and 10000 r g suffices 
to cover the observational data. The simulated core red- 
shift values follow a power law, z = pr s , for fixed incli- 
nations and are plotted in Fig. [7| Observational data for 
Mrk 110 by K03 are overplotted as boxes. The horizontal 
error bars are due to the uncertainty in time lag measure- 
ments and the vertical error bars involve the uncertainties 
in the differential redshifts of the rms line centers. The 
labels of the axes are adapted to theoretical considera- 
tions and display redshift z as a function of distance to 
the black hole in units of gravitational radii. To rescale 
the x-axis the best-fit value for the black hole mass 2 ob- 
tained by K03 is assumed, M ~ 1.4 x 10 8 M . Tab. 
shows the fitting results for power laws at each inclina- 
tion angle. At any inclination the power law exhibits the 
same slope, s num w —1.002 ± 0.005, which is the average 
of the five inclination angles assumed here. This is a di- 
rect consequence of the Schwarzschild factor, z(r) oc r" 1 
i.e. Sthco = — 1- Additionally, the power laws shift toward 
lower gravitational redshifts as the inclination angle in- 

2 In general, it is recommended to treat M as free parameter, 
too. 
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Table 2. Power law fit parameters p and s for inclinations 
1° to 40° 





P 


error Ap 


s 


error As 


i 


1.566 


± 0.007 


-1.001 


± 0.001 


10 


1.428 


± 0.008 


-0.999 


± 0.001 


20 


1.187 


± 0.003 


-0.996 


± 0.001 


30 


0.914 


± 0.013 


-1.007 


± 0.003 


40 


0.484 


± 0.013 


-1.006 


± 0.005 



KBHRTi =01 
KBHRT i = 10 
KBHRTi =20 
KBHRTi =30 
KBHRTi = 40 
Kollatschny 2003 I 
best fit 0.88/r 




creases. This is due to projection effects - included in the 
projection parameter p. The inclination dependence of the 
redshift can be approximated by z oc cos(i) i.e. p is related 
to the cosine of i. In other words, the higher the inclina- 
tion, the more blueshift there will be (see Fig.[3J), and this 
Doppler blueshift counteracts the redshift. Interestingly, 
this behaviour could be exploited to determine the inclina- 
tion angle of the inner disk from observed gravitationally 
redshifted features. 

The best-fitting power law is shown as a thick solid 
line in Fig. Here, the slope of s = —1 has been fixed 
and the projection parameter p fitted to give PMrkiio — 
0.886 ±0.034. Its relation to a specific inclination i can be 
extracted from Fig. [3] which shows the cosine behaviour of 
p for the simulated sets i S {1, 10, 20, 30, 40°}. The cosine 
fit to ray tracing data yields p(i) ~ 4.63 cos(i) — 3.13 (error 
~ 10%; note that this is only valid for i < 40°) within the 
radial range between 100 and 10000 r g . Taking the best 
fit value 0.886 one reads in Fig. |H] at the cross that the 
inner disk of Mrk 110 is inclined to i ~ 30°. This result is 
consistent with that of K03 (i ~ 21 ± 5°). 

4.3. Static vs. stationary emitter velocity field 

The redshifts calculated above were based on a simple, 
but dynamical BLR model. In this section we clarify un- 
der which conditions the Schwarzschild factor can be used 
to estimate gravitational redshifts. Usually, a velocity field 
of the BLR emitters has to be assumed. The prominent 
Schwarzschild factor can only be applied to quantify the 
redshift in case of static emitters (and static black holes). 
Of course, real emitters such as the BLR are dynamical 
e.g. in stationary motion. In this case, the ray tracing 
technique is a convenient method for redshift computa- 
tions. This approach is also justified by the fact that it 
takes into account that real emitters are extended. The 
relativistic <?-factor (consult M04 for details) for static 
emitters simplifies significantly because the Keplerian an- 
gular velocity satisfies Six = for radii larger than the 
radius of marginal stability, r > r ms . The emitter velocity 
field in Bardeen Observer's frame is assumed to be purely 

rotational i.e. v^' = with i = r, 6. We are thus left 

r /2 

with a Lorentz factor 7 = (l — V(^\V^') . Considering 
fi = Q K = in v m = Co (^ L ) the Lorentz factor simply 
goes to unity, 7 — > 1. It should be kept in mind that this 



io-= 1 ■ ■ ■ — ■ ' ■ — 1 

100 1000 10000 

distance [rj 

Fig. 7. Radial dependence of the line core redshift: data 
for the NLS-1 Mrk 110 taken from K03 (boxes) are com- 
pared to the line core redshifts as computed from Kerr 
ray tracing simulations with different inclination angles. 
Redshift as a function of distance scales with a power law 
with identical slope s = — 1 for all inclinations but with 
different projection parameter p that determines the ver- 
tical shift of the power law. Best fit for Mrk 110 (solid 
thick) yields p ~ 0.886 that can be used to determine the 
inclination of the inner disk. 
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i = 29.7 for Mrk 110 





































10 20 30 40 

inclination angle [degree] 

Fig. 8. Variation of the projection parameter p with incli- 
nation i as computed from Kerr ray tracing. The numeri- 
cal data points can be approximated by a cosine function 
(solid), p(i) ~ 4.63 cos(i) — 3.13 within 10 % precision, p 
measurements from observations give inner disk inclina- 
tions. The cross marks the result for Mrk 110 based on 
K03 data: i ~ (29.7 ± 1.2)°. 



regime is not applicable at r < r ms when frame-dragging 
u> becomes important and f2 = is replaced by £1 — ► oj. 
Altogether, the generally complicated <?-factor is identical 
to the lapse function a for static emitters in the r > r ms 
regime. The quantity a is evaluated here for static black 



8 



A. Miiller and M. Wold: Signatures of Gravitational Redshift 



holes (as already shown in K03) and for rotating black 
holes. Restriction to the equatorial plane, 8 — ir/2, yields 
a simple expression for the lapse function (see e.g. M04, 
Eq. 3): 



r(r 2 ~2Mr + a 2 ) 
a ~ V r 3 + a 2 (r + 2M) ' 

The lapse function reduces 
Schwarzschild factor for a = 



(4) 



to the well-known 



a L=o= 



2M 



(5) 



4.4. 6/ac/c /70/e rotation 



A rotating SMBH can not be excluded a priori in the 
case of Mrk 11 0. On the cont rary, theories of black hole 
growth, see e.g. lShapirol (|20f)5l > . strongly suggest fast spin- 
ning SMBHs in the local universe with either a/M ~ 0.95 
(MHD disk) or even a/M ~ 1 (standard thin gas disk). 
However, it is also known that the rotation of space-time 
decays very rapidly in the outer regions of the black hole 
gravitational potential. The frame-dragging frequency de- 
creases according to u> oc r~ 3 . This is documented in the 
final Fig. which summarizes the results from Fig. [7] and 
IHI and also extends beyond the region explored by K03. 
The redshifts corresponding to g-factors of g i — ► a are 
plotted for a non-rotating and a fast rotating black hole 
with a = (solid curve) and a/M = 0.998 (dotted curve), 
respectively. The curves are compared to the observations 
of Mrk 110 (boxes) and the best-fitting i — 30° Kerr ray 
tracing simulation (filled circles). It is evident from the 
figure that the rotation of space-time is important only 
at small radii, r < 4 r e . Therefore, black hole rotation 
can only be probed with spectral features originating in 
regions very close to the black hole, like the X-ray fluores- 
cence lines of iron. Optical emission lines are not suited for 
probing black hole rotation in that manner - at least not 
for AGN. The reason for the small offset between the ray 
tracing results shown by the filled circles in Fig. and the 
Schwarzschild/Kerr lapse functions plotted as solid and 
dotted lines is the different velocity field of the emitters 
and the transverse Doppler effect. The ray traced emitters 
are assumed to rotate stationarily, whereas the lapse func- 
tions only coincide with ray traced (/-values when static 
emitters are assumed. Also the ray traced redshift values 
are based on g core computed by averaging over the whole 
line (Eq. |2J) , and therefore are less sharp than those given 
by the lapse function. 

Multi-wavelength observations such as the ongoing 
COSMOS project 3 may help filling the gap between 2 and 
200 r g in Fig. El If several gravitationally redshifted spec- 
tral lines can be identified in a source, the central mass 
can be determined with high accuracy. 



KBHRT i = 30 
Schwarzschild lapse 
Kerr lapse 
Kollatschny 2003 




http : //www . astro . caltech . edu/ ~cosmos/ 



1 10 100 1000 10000 

distance [r g ] 

Fig. 9. Synoptical plot with optical K03 data for Mrk 110 
(boxes) and best fitting Kerr ray tracing simulation with 
i ~ 30° (filled circles) as well as lapse functions for a static 
(Schwarzschild, solid) and rotating (Kerr, dotted) black 
hole. An essential statement illustrated here is that optical 
BLR lines can not probe black hole rotation due to their 
huge distance. Generally, multi-wavelength observations 
are recommended to fill the gap at smaller radii. 



4.5. Assumptions 

BLR emission lines are influenced by both gravitational 
redshift and Doppler effects. Therefore, the fitting pro- 
cedure involves two parameters, black hole mass M and 
inclination angle of the emitter i. The inclination can be 
determined if the black hole mass is known from other 
methods. In case of Mrk 110, the mass of the SMBH was 
inferred from reverberation mapping as well as from grav- 
itational redshift in the literature. Other techniques are 
measurements involving stellar and gas dynamics, the M— 
a relation and Masers. To be able to determine the inner 
inclination, we used the best fit black hole mass computed 
from gravitational redshift found in K03. The error boxes 
of K03 can be fitted quite well with the ray tracing results 
(3.9 % error). Hence, a model with stationary rotating 
rings distributed in a radial range between 240 and 4700 
r g and inclined to i ~ 30° can explain the observed gravi- 
tational redshift of optical lines in Mrk 110. It is surpris- 
ing that a simple Keplerian rotating model can describe 
the BLR structure so well . Similar results have been ob- 
tained e.g. for NGC 5548 llPeterson fc Wandellll999h and 
other AGN llPeterson fc WandefeOOOHOnken fc Peterson! 
2002; Chen et al. 1989). In this work we have made the as- 
sumptions of a flat and Keplerian rotating BLR structure. 
However, nature is supposed to be more complicated al- 
lowing particularly for BLR wind components and o pacity 
effects IjPopovic et al ll995UChiang fc Murravll 99~1 . This 
implies a modification of the solver including radial and 
poloidal motion as well as photon propagation through 
a dense volume. Furthermore, disk warping has not been 
considered and this may significantly change the results. 
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For warped disks the analysis has to be extended by using 
suitable generalize d ray tracing codes e.g. like the one by 
ICadez et all Ir2003h . 

5. Conclusions 

Line cores at distances from 2 to 100000 r g from 
a rotating black hole have been analysed using rel- 
ativistic ray tracing simulations in the Kerr geome- 
try. The line cores are gravitationally redshifted by 
z coro ~ 10~ 5 , 10~ 4 , 10~ 3 , 10~ 2 , lO" 1 10° at distances of 
100000, 10000, 1000, 100, 10, 2 r g from the black hole, re- 
spectively. This z ps dista „ cc j r j behaviour at large radii 
is a straightforward consequence of the Schwarzschild fac- 
tor. Lines characterised by a core energy <7 coro confirm this 
scaling behaviour. 

Gravitational redshift occurs in two modes. One 
regime starts at larger distances from the black hole and 
shifts only the line as a total feature while conserving its 
intrinsic shape; the amount of redshift can be looked up 
in Tab. ^ The other regime, which is the strong gravity 
regime, dominates at r < 10 r g or at r < 10 13 cm for 
a ten million solar mass black hole. Relativistic emission 
lines originating in this region are strongly deformed and 
suppressed and differ substantially from the corresponding 
line profiles in the emitter's rest frame. 

The present work has demonstrated that the onset of 
GR becomes important at distances smaller than 75000 r g 
if assuming an optical spectral resolution of 0.1 A. It is 
stressed that this critical radius depends on the astro- 
nomical resolving power and lies farther out if the resolu- 
tion is higher. However, even with high-resolution spec- 
troscopy it remains a challenge to probe gravitationally 
redshifted spectral lines due to the fact that competing 
effects and more complex physics are likely to be involved: 
a flat Keplerian BLR model may be too simple to model 
the complex BLR velocity field and both spherical BLR 
structure and a wind component may play a crucial role. 
In addition, narrow-line components and e.g. contamina- 
tion of the H/3 line by Fell may also complicate both the 
analysis and the interpretation. Nevertheless, the theory 
of General Relativity predicts that the gravitational red- 
shift is present, and a valuable ansatz for probing it is to 
search for systematic shifts with varying distance, as done 
by K03. It is suggested in this work to supplement this 
by multi-wavelength observations that should all point 
towards the same central mass and inner inclination. 

We confirm the analysis by K03 of the NLS-1 galaxy 
Mrk 110 here, using a more general treatment with sta- 
tionary (/-factors that are in concordance with observa- 
tional data. Ray tracing simulations in the Kerr geometry 
support an inclination of i = 30° for the inner disk of 
Mrk 110. Reversely, if the inclination is known, this can 
be exploited to determine the black hole mass from the fit- 
ting procedure outlined here. Whether fitting i or M (or 
both) - such techniques may help explore AGN unification 
schemes: multi- wavelength studies allow for studying the 
inclination deep into the AGN and to probe orientation 



and luminosity-dependence of AGN types. Furthermore, 
we show that broad optical lines can not serve as a probe 
of black hole spin because frame-dragging effects only oc- 
cur very close to the black hole. In this region, only hot 
emission lines (such as the Fe Ka line in X-rays) or other 
relativistic spectral features indicate black hole spin. The 
analysis presented here is not only valid for supermassive 
black holes but also for stellar-mass black holes in X-ray 
binaries or intermediate-mass black holes that may be 
found in ultra-luminous X-ray sources or globular clus- 
ters. 
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